Hamiltonian decomposition for bulk and surface states 
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We demonstrate that a tight-binding Hamiltonian with nearest- and next-nearest-neighbor hop- 
ping integrals can be decomposed into bulk and boundary parts in a general lattice system. The 
Hamiltonian decomposition reveals that next nearest-neighbor hopping causes sizable changes in 
the energy spectrum of surface states even if the correction to the energy spectrum of bulk states 
is negligible. By applying the Hamiltonian decomposition to edge states in graphene systems, we 
show that the next nearest-neighbor hopping stabilizes the edge states. 



The energy band structure is of central importance in 
understanding the electronic properties of material. A 
tight-binding (TB) model is a versatile approach to study 
the electronic, magnetic and transport properties of solid 
since TB model describes the qualitative features of the 
energy band structure. [l| The long-range hopping terms 
such as next nearest-neighbor (nnn) hopping are often 
added to the TB model with nearest-neighbor (nn) hop- 
ping to improve the energy band structure. In many 
cases, the nnn correction to the energy band structure is 
not a matter of particular importance. In this Letter, we 
show that the nnn hopping can change appreciably the 
energy spectrum of surface states which appear near the 
boundary of a system even when the correction to the 
energy spectrum of bulk states is negligible. We explain 
this by decomposing the nnn TB Hamiltonian into two 
parts; bulk and boundary parts. This Hamiltonian de- 
composition is essential to understanding the stability of 
surface states. 

We use graphene systems to demonstrate the Hamil- 
tonian decomposition for the following reasons. (1) 
Graphene is known to have both bulk and surface states. 
The bulk states exhibit a "relativistic" energy band struc- 
ture called the Dirac cone, 0, Q and the surface states 
called the edge states appear near the zigzag edge. |4j 
(2) The edge states have been observed by several ex- 
perimental groups, [1, @, 0, El and a theoretical under- 
standing of the experimental results is called for. In fact, 
the nnn hopping is important to explain the experimen- 
tal results. [9( (3) The edge states have a large density 
of states (DOS) near the Fermi energy, which is respon- 
sible for the Fermi instabilities. Since the DOS depends 
on the the energy spectrum or bandwidth of the edge 
states, the nnn hopping is important for the appearance 
of many-body effects of the edge states. [l(J U| Thus, 
graphene is a good testing system for the Hamiltonian 
decomposition not only from a theoretical but also from 
an experimental point of view. 

We study the nn (nnn) TB Hamiltonian, H nn (H nnn ), 



which is defined as 



H, 



"7n 



= c\ [7i nnn ] ij Cj , 



(1) 



where Cj (cj) is the annihilation (creation) operator of 
an electron at i-th site, 70 (7,1) is the nn (nnn) hopping 
integral, and the matrix element [7i n n]y ([Hnnn]ij) is 1 
when i-th site is a (next) nn site of j-th site and is zero 
otherwise. In the following, we will show that H nnn can 
be decomposed into bulk and boundary (edge) parts for 
a graphene with zigzag edge. 




FIG. 1: (a) An electron at the central site, \^o) = c o|0), is 
transfered to three nn sites by H un /{— 70). The resultant state 
is |*o) = Kigali [Wnn]ioc||0). (b) H nn transfers the electron to 
the nnn sites. At the same time, H nn returns the electron to 
the original site. The matrix element of {H. nn ]ij that returns 
the electron to the original site is given by 3 because there are 
three nn sites around the central site, (c) The matrix element 
for the return process [W„ n ]ij depends on whether the site is 
a bulk site ([TC nn ]n = 3) or a zigzag edge site ([H nn ]jj = 2). 

Suppose that we put an electron on the central site 
denoted by the empty circle in Fig. QJa). The initial 
state is labeled as |\& ) = Cq 1 0) . We operate on |^ ) 
with H nn / (—70), then H nn transfers the electron to three 
nn sites denoted by the solid circles in Fig. [TJa). This 
state is written as \^' ) — J2iea.iiP^nn]io c l\0) ■ The num- 
bers associated with the solid circles in Fig.QJa) indicate 
the matrix element of Ti nn . The successive operation of 
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£W(-7o) on |*o) gives 

K) = E [Wnnk[W»n]«)ct|0). (2) 
i, j Gall 

Starting from the initial site, the electron reaches the 
nnn sites as shown in Fig.[ljb). Thus, the two successive 
nn hopping processes relate to the nnn hopping process. 
This indicates that H nn includes H nnn . However H nn 
and 7~C nnn are not identical because in H nn there is a 
diagonal matrix element that returns the electron to the 
original site, that is, |\l/g) of Eq. {2| contains the term 
with j = 0. Since there are three nn sites around the 
original site, the amplitude of this return process is 3 as 
shown in Fig. [T^b). H nn and Hnnn become identical if 
we subtract the corresponding diagonal matrix clement 
from H nn . This matrix is proportional to the unit matrix, 
[T\ij = Sij. 

For a periodic system, since the number of bonds of 
every site is three, we have H nnn = H nn — 3/. Putting 
this into Eq. |T]), we see that H nnn can be rewritten as 

H nnn = -J n c ll W nn — 37]yCj. (3) 

ij'Gall 

The matrix H n n can be diagonalized by a unitary matrix 
as [UH nn U^] pq = E p /(~^ )S pq: where E p is the energy 
eigenvalue of an eigenstate \E V ). Then, from Eq. (|3|) we 
see that H nn and H nn n can be diagonalized simultane- 
ously by the basis of \E p ), and the energy eigenvalue of 
the total Hamiltonian, H nn + H nnn , is given by 

for \E p ). In Eq. ©, H nnn contains the on-site potential 
part, 37 n X^ieaii C I C »- This on-site potential can be ig- 
nored since it changes only the origin of the energy band 
structure, as shown by 3j n in Eq. (j4|). 2] This statement 
is correct for a periodic system without boundary, but is 
not approved for a system with boundary. It is because 
of that the number of bonds of the edge sites is different 
from that of a bulk site and the corresponding on-site 
potentials at the edge sites are different from those at 
the bulk sites. 

To show this explicitly, we put an electron on the zigzag 
edge site labeled as 2 in Fig.QJc). The electron is trans- 
fered to the nnn sites by H^ n . For this time, however, the 
matrix element of the on-site potential term that we need 
to subtract from H nn in order to get Hnnn is 2 because 
the number of bonds is 2 for the edge site. It is different 
from 3 for a non-edge (bulk) site. Thus, we obtain the 
formula for H an n as 

Hnnn = -Jn ^ C H[^lnh ~ 9i > ' i\<r (5) 
ij'Gall 



where is the number of bonds of i-th site. Since we 
can shift the origin of the energy by 37„ without a loss 
of generality, Eq. ((5|) can be written as 

Hnnn = ~ln ^ 4W-nn]ij c i + In ^ _ 3 )^' ( 6 ) 
zjGall iGall 

where hi = ejej is the number operator of i-th site. The 
first term in Eq. ([6]) (or the second term of Eq. ((4])) shows 
that Hnnn breaks the particle-hole symmetry of H nn be- 
cause Hnnn contains the square of Hnn- The second term 
in Eq. (JSJ) represents quantum well potentials at the edge 
sites with potential depth of —"f n because gt = 2 for a 
zigzag edge site. The quantum well potential appears 
only at an edge site whose number of bonds is differ- 
ent from those of a bulk site. Therefore, a surface state 
appearing near the boundary is strongly affected by this 
potential. For bulk states, only the first term of the right- 
hand side of Eq. ([6]) is important. In fact, if a system has 
no boundary (if a system is periodic) , the second term of 
Eq. |6]) disappears and the nnn Hamiltonian is given only 
by the first term. Therefore the first term can be consid- 
ered as the bulk part and the second term is as the edge 
part of the nnn Hamiltonian, i.e., H nn n — H^^ + H^f^, 
where 

ffbulk = i 

ijGall 

ieall 

The edge state is labeled by the wavevector along the 
zigzag edge, k, as \E(k)). Here E(k) denotes the en- 
ergy eigenvalue of H nD - Since E(k) of the edge state 
is very close to zero, [J] the energy shift due to H^ n k , 
— {lnha)E{k) 2 1 is negligible. Thus, the energy correc- 
tion to the edge state arises from H nn n as 

AE(k) = -7„ E (E(k)\hi\E(k)). (8) 

iGcdgc 

It is only the density at the edge sites that determines 
AE(k). The energy bandwidth (W) of the edge states 
can be calculated in the following way. As k approaches 
to the Fermi point, i.e., ka — * 27r/3 or Att/3 (a is lat- 
tice constant), the edge state changes into a bulk state 
since the localization length £(fc) — > oo. 0, [12] Then 
E Jeo dgo(£( fc )N£( fc )> is negligible and A£(2^/3a) = 
AE(An/3a) = 0. On the other hand, the k = n/a 
state is the most localized state satisfying £(ir/a) = 0. 
For this state, we have J2iecd S c(E(k)\ni\E(k)) = 1 and 
AE(tt/o) = — 7„. It shows that the k — ir/a state (the 
Fermi point) is located at the bottom (top) of the en- 
ergy band and W — j n for the edge state. This re- 
sult of W = 7„ reproduces the energy bandwidth that 
is numerically calculated for the zigzag edge shown in 
Fig. [U [9( The quantum well potential of ff^l lowers 
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the edge state's energy if j n is a positive value. We adopt 
7„ w 0.3 eV. This value is obtained by a first-principles 
calculation with the local density approximation. [13j 




FIG. 2: The energy band structure of graphene system with 
the zigzag edge. This plot is obtained by diagonalizing H nn + 
H nnn — 37„ numerically. We adopt 70 = 3.0 eV and y n = 0.3 
eV. The horizontal axis is a wavevector along the zigzag edge 
(k) multiplied by the lattice constant (a). 

It should be mentioned that Eqs. ([5]) and © include 
nnn hopping between sites which are connected by two 
successive nn hopping processes. Two successive nn hop- 
ping processes do not include disconnected nnn hopping 
process. The disconnected nnn process is relevant to the 
Klein edge as shown in Fig. [3Ja). [3] Because H^ n can 
not transfer an electron at the edge site (i) to the nnn 
edge site (j or j') in Fig. 02(a), we have to add H„ nn , 
which represents the nnn hopping between disconnected 
sites, to the matrix of the right-hand side of Eq. ([5]) in 
order to get a complete Ti. nnn . Thus we have 

ffnnn = ~7n E 4{[^nJy ~ 9i[I\ij + [^nL]«} c i' 

(9) 

Since [7i„ nn ]ij * s n °t zero only when i and j are both the 
edge sites, the disconnected nnn Hamiltonian is written 
as 

^111111 = ~ln ^ t c i [W nm ] i j Cj 
nnn 

= - 7n Yl c \ c r ( 10 ) 

i,jEcdgc 

H^nn can be classified into the edge part of the nnn 
Hamiltonian since H^ n is given by the creation and an- 
nihilation operators at the edge sites. If we represent the 
wavefunction using the density n and phase 9 as 

i*>=e wmi^ (n) 

igall 

then we have 

nnn 

<*|fl£n|*> = -7» E y/m^- *- (12) 

ij'Ecdgc 



This result shows that not only the density (y/niUj) but 
also the relative phase (0* — 0j) of the localized wavefunc- 
tion is important for the energy shift. This is contrasted 
to the fact that the quantum well potential couples only 
to the density of a quantum state. If there is a lattice 
periodicity along the edge, we can set 0j = (ka)i and 
rij = nj. Then Eq. (|12p becomes 

<*|ff£nl*)=-27nCO S (H Ui. (13) 

iGcdgc 




FIG. 3: (a) The lattice structure of the Klein edge. The nnn 
hopping between nnn edge sites (i and j or i and j') is not 
represented by the double of the nn hopping, (b) The energy 
band structure of graphene system with the Klein edge. 

Since g L — 1 for the Klein edge sites, the edge part of 
the nnn Hamiltonian is written as 

nnn 

iGcdgc i,j Gcdgc 

Using Eq. (H3D, we get AE(k) (= {E(k)\H°*f\E(k))) for 
the Klein edge as 

AE(k) = -2 7 „(l + cos(fca)) E n *( fc )' ( 15 ) 

iGcdgc 

where 7ij(fe) = (£ , (fc)|ni|£ , (fc)). The energy bandwidth 
for the Klein edge states is calculated as follows. Near 
the Klein edges, the edge states appear for < k < 2ir/3a 
and 47r/3a < k < 2tt/o, (see Fig. EJb)). It can be shown 
that the wavefunction of most localized state is given 
by k — state, and sum of the densities at the Klein 
edge sites is given by Z^gedge = 3 / 4 - Ell Ther h by 
putting k = into Eq. (fl5|. we have AE(Q) — — 37„. 
Thus, W for the Klein edge is 37„ (=0.9eV). This result 
also reproduces the energy bandwidth that is numerically 
calculated for the Klein edge shown in Fig.[3Kb). The half 
of W comes from the quantum well potential and the rest 
half of W is due to the disconnected nnn hopping process. 

Here let us summarize the formula for the nnn Hamil- 
tonian: -ffnnn can be decomposed into bulk and edge 
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parts as H nnn = + with 



edge states. A theoretical calculation of the self-energy 



rrbulk _ _ _tr<lj2 1 

n niin — In / , c i L rl nnJ*j l 
i, j Gall 



iGcdgc i, j £edge 



(16) 



g is the number of bonds of a bulk site. The first term of 
the right-hand side of H^f° in Eq. (|T5]) represents quan- 
tum well potentials at the edge sites. It is only the num- 
ber of bonds at the edge site which determines the depth 
of the quantum well potential. As a result, the quan- 
tum well potential appears at the edge sites, regardless 
of the edge shape. For example, in the case of graphcnc, 
9i 7^ <?(= 3) not only for the zigzag (or Klein) edge sites 
but also for the armchair edge sites. Thus, for a finite 
system of graphcnc shown in Fig. |4j the quantum well 
potentials of —j n are denoted by the solid circles. The 
quantum well potential at the armchair edge may be not 
of importance as the zigzag edge because the edge state 
of a graphene is absent from the armchair edge. 0, @, @] 
It is also interesting to note that g% ^ g for the three 
sites around a lattice vacancy. In addition to the quan- 
tum well potentials at the three sites, the disconnected 
nnn hopping appears between them. 



for the edge states is given in Refs. [15| and [11 . 

We note that Eq. (fl"6)) is not restricted to graphene sys- 



tems but is applicable to other two-dimensional systems 
like the square lattice, and three-dimensional systems. 
The nnn Hamiltonian stabilizes surface states in a general 
system through the quantum well potentials. The energy 
dispersion relation of surface states is observed below the 
Fermi level by high-resolution photo-emission studies of 
the (111) surfaces of copper, silver and gold. [l(| We spec- 
ulate that the observed stability of the surface states is 
due to the edge part of the nnn Hamiltonian. 

In conclusion, H nnn can be decomposed into the bulk 
and edge parts as shown in Eq. (fTS)) . If the energy spec- 
trum of H nn is symmetric with respect to E — 0, then 
Hj* breaks this symmetry. If H nn has a localized en- 
ergy eigenstate near the edge of a system, then -ff„ n f c is 
relevant to shift the energy eigenvalue through the quan- 
tum well potential and the disconnected nnn edge Hamil- 
tonian. The quantum well potential couples only to the 
density, whereas the disconnected nnn edge Hamiltonian 
couples to the phase of the localized wavefunction. Al- 
though the Hamiltonian decomposition is proved for two- 
dimensional graphene systems, it works for other systems 
like a one-dimensional chain of atoms and a three dimen- 
sional lattice system as well. 
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FIG. 4: A graphene system with a boundary and a lattice 
vacancy. H^tf-? gives an on-site potential energy shift of — j n 
at carbon atoms denoted by solid circles. There is a discon- 
nected nnn Hamiltonian at the dotted lines around the lattice 
vacancy. 

We remark on other effects that can modify the energy 
spectrum of the edge states. First, it is naively expected 
that the orbital energy at an edge carbon atom is differ- 
ent from that at a bulk atom when a functional group at- 
taches to the edge atom. The attachment of a functional 
group gives rise to an additional change of the energy 
bandwidth of the edge states. However, this energy shift 
may be positive or negative value depending on the type 
of a functional group. This can be distinguished from 
the energy shift due to the nnn Hamiltonian because it 
is always a negative value. Second, W can be modified 
by the electron-electron or electron-phonon interactions 
because they give rise to a self-energy correction to the 
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